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1. Introduction 

Polarization is a fundamental property of electromagnetic fields. Its effects are partic- 
ularly important when the field's 3-dimensional, vectorial character plays a role, such 
as radiation from sources, the optical near-field, tightly- fo cussed beams, and scatte- 
ring. In these situations, the conventional theory of partial paraxial polarization fails. 



There has been recent interest in extending the theory of partial polarization to these 
nonparaxial situations, particularly the notion of degree of polarization (e.g. Refs. 
1-7 and references therein) . 

The transverse polarization properties of optical beams are well understood. ^'^ 
The theory is based on the hermitian 2x2 polarization matrix p^-, defined 

_ {E*E,) 

with averaging over some ensemble of transverse field vectors E — {Ei,E2). This 
matrix has a unique linear decomposition: p2 = -P2Ppure + (1 ~ -^2)P2,un) where Pp^re is 
the polarization matrix (single-point coherence matrix) of a pure polarization state, 
P2,un is the trace-normalized 2x2 identity matrix (completely unpolarized) , and 
< P2 < 1- -P2 is the modulus of the difference between the eigenvalues of P2, so 

2 _ 2trp|- (trp2)' 

When P2 is decomposed via the basis of Pauli matrices, the coefficients are the Stokes 
parameters 5*1, 5*2, 5*3 (5*0 = 1 since tr/Og is normalized). The sum of squares + 
82+3^ = P2] on propagation through rotators and retarders (represented by unitary 
transformations on pg); the Stokes parameters may change, but P2 does not. 

The theory of 3-dimensional partially coherent fields is less well developed. Even 
in fully polarized fields, the Stokes parameter description fails since the plane of the 
polarization ellipses varies with position. Furthermore, in partially polarized fields, 
polarization states with ellipses in different planes may be incoherently mixed. Sev- 
eral measures of the 3-dimensional degree of polarization have been proposed and 



discussed;^' ^"^ these measures, discussed below, have the property of unitary invari- 
ance: they only depend on the eigenvalues of the 3-dimensional polarization matrix 
P3, by analogy with unitary invariance of P2- In terms of a measure on the eigenval- 
ues, P2 is clearly unique (up to monotonic transformations, such as squaring), since 
after trace normalization, the two eigenvalues have only one freedom (their differ- 
ence). However, for 3 x 3 polarization matrices, there are two such freedoms, and 
there is no unique mathematical definition for the degree of polarization; it is there- 
fore appropriate to examine measures that emerge naturally out of physical processes, 
even if desirable mathematical requirements have to be relaxed. Clearly, one such re- 
quirement that cannot be satisfied in three dimensions is the decomposition of the 
polarization matrix into purely polarized and purely unpolarized parts;^ the inter- 
mediate eigenvalue must also play a role in any matrix decomposition (the largest 
and smallest eigenvalues corresponding to the completely polarized and unpolarized 
parts) . 

This is the optical analog of a well-known problem in quantum mechanics, namely 
the non-uniqueness of the decomposition of a general density matrix. Thus, the 
definition of the degree of polarization is chosen here to correspond to that of a purity 
measure in quantum mechanics: a purely polarized state (with eigenvalues 1,0,0) 
always has measure 1, and the completely unpolarized state (1/3 times the identity 
matrix), isotropic in all directions, has measure 0; any other state of polarization has 
a value in between. For P2 clearly satisfies this, and is effectively unique. However, 
under this definition, the quantity proposed in Ref. 4 (the difference between the two 
larger eigenvalues of pg) - motivated by the eigenvalue decomposition of polarization 
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- is not a purity measure, as anisotropic states of 3- dimensional polarization, which 
have directional information, are counted as completely unpolarized (for instance, a 
completely unpolarized paraxial beam has directional information (its propagation 
direction), but its eigenvalues are 1, 1, 0). In quantum mechanics and paraxial optics, 
such unitary invariance is a physical requirement, but there seems to be no strong 
argument for this in the nonparaxial point discussed towards the end of this 

article. 

Here, I describe an alternative measure for the degree of polarization in three 
dimensions, based on isotropic Rayleigh scattering, which is not a unitary invariant 
of P3. This is defined and discussed in the following section, then in Section [3], is 
compared with various other 3-dimensional polarization measures for various specific 
examples of P3. The final section consists of a discussion of various issues associated 
with the Rayleigh-defined degree of polarization and non-unitarity. 

2. The Rayleigh-defined degree of polarization 

Rayleigh scattering is a fundamental 3-dimensional optical phenomenon, and is of 
primary importance in scattering theory. ^^"^^ At an isotropic scatterer, of dimension 
much smaller than the optical wavelength, the possibly incoherent electromagnetic 
field is represented by density matrix P3 (defined analogously to p2 in Eq. ([1])). Par- 
tially polarized rays are scattered in all directions 6', 0; the properties of the scattered 
ray in 9, are determined by the polarization matrix^^'^^ 




(3) 
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where a is the total scattering cross-section, and p_l(6', 0) is a projection matrix, 



projecting into the plane perpendicular to the unit vector u{9, (p) in the direction 
9, (p. Although p2{9,(j)) is 3 X 3, it is null in (appears transverse to) the direction 
u{9,(j)) by definition. p2{6,(j)) is not trace-normalized, its trace being the intensity 
1(9, (h) scattered in the 9, (/)-direction: 



I{9, 



tr p2( 



(4) 



(this is the differential scattering cross section^^'^^). The scattering mechanism is 
here idealized such that the scatterer experiences no recoil (there is no center-of-mass 
motion in the scattering). 

The degree of polarization of the scattered light depends on the scattering direc- 
tion, but is otherwise similar to Eq. ([2]), 

2tipli9,<P)-[tip2i9,<P)]' 



2 ' 



(5) 
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where the denominator equals I{9,(j)). The Rayleigh- defined degree of polarization 
-pRay is taken to be the average over all directions of P{9,(f)), weighted with respect 
to I{9,(f)) : 



p ^ lo d<l>j;d9 sin gP2(g, 0)^0) 
''''' J^''dcj)J^d9sm9Ii9,(P) 

where the denominator is the total scattering cross-section a. 

This process can be visualized using the geometric interpretation of the 3- dimen- 
sional polarization matrix.'^ Choosing appropriate cartesian axes, p^ can be written 

Ml -iN3 iN2 
P3 = iNs Ms -iNi , (7) 



V 



-iN2 iNi Ma 



where Mj.Nj are real, and 1 > Mi > M2 > M3 > 0. The real part, M = 
diag{Mi, M2, M3}, is interpreted as the moment of inertia ellipsoid of the polarization 
ellipse ensemble. The imaginary part is an axial vector, N = [Ni, N2, N^), the aver- 
age angular momentum vector of the ensemble. A completely unpolarized ensemble 
(representing, for instance, black-body radiation) is therefore represented by a sphere, 
and a pure polarization state by an ellipse with orthogonal vector (whose length is 
the ellipse area). For a pure polarization state, the 1 (2)-axis is the major (minor) 
axis of the polarization ellipse, and the 3-axis is its angular momentum direction: the 
fiat ellipsoid (M3 = 0) has axes proportional to the polarization ellipse axes squared 
{Mj = \E^\'^/\E\'^,j = 1,2) and \N\ = {N^l = \E* x E\/2 = ^MX. The ellipsoid 
and vector can change under unitary transformation. 

The ellipsoid M and vector N generalize features of the real part Re p2 de- 
fines an eUipse, and the imaginary part, proportional to the Stokes parameter S'3, 
is the overall angular momentum perpendicular to the transverse plane. The polar- 
ization state P2{0,(f)) of a scattered ray is geometrically determined by projection; 
the transverse ellipse of its real part is the projection of the ellipsoid in the plane 
perpendicular to 6,(f), and its imaginary part is the projection of N into u{6,(f)), as 
depicted in Fig. [1] projecting in the coordinate axis directions. Paay is the average of 
the degree of polarization over all these projections. 

Although -Paay depends only on the polarization matrix pg, it cannot be simply 
expressed in terms of the matrix elements. Nevertheless, it is a purity measure on the 
density matrix pg, as defined above. If pg represents a pure state, its eigenvalues are 
1,0,0, and so P2{6,(f)) (the difference divided by the sum of eigenvalues of p2{6,<f))) 
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Fig. 1. Geometric representation of polarization matrix P3 as ellipsoid and 
vector. The ellipses and vectors for p2{6,(f)) in three orthogonal projection 
directions are also shown. 

is 1 for almost all directions, thus Paay = 1- Conversely, Paay = 1 in Eq. ([6]) only if 
P {$,(/)) = 1 for almost all directions, which only occurs if one eigenvalue of p2{0,(j)) 
is zero, implying pg is pure. Similarly, -pRay = if and only if pg is the completely 
unpolarized 3-dimensional matrix 1/3. 

In this argument, 'almost all' assumes its technical sense, that is, for all values 
contributing to the integrals in ([6]). If P3 represents a pure state of linear polariza- 
tion diag{l, 0, 0}, the ray scattered in the 1-direction has zero intensity, but such an 
isolated direction does not affect -pRay 

3. Comparison with other measures of degree of polarization, and values 
for sample polarization matrices 

In this section, specific examples of 3-dimensional ensembles represented by particular 
polarization matrices pg are considered, demonstrating explicitly that Paay is not a 
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unitary invariant; progress by means of example is the only way understanding of PRay, 
in the absence of a general closed analytic form. Before discussing these examples, I 
will describe some other measures which are unitary invariant, that is, they depend 
solely on the eigenvalues 1 > Ai > A2 > A3 > 0. The measure of Ref. 4, defined 
Ai — A2, is not a purity measure of polarization under the present definition; however, 
a slightly adjusted definition is: 

-Plin = Ai — A3, (8) 

which is obviously linear in the eigenvalues. A second measure, introduced in Ref. 2, 
depends on a quadratic combination of the eigenvalues, defined 

^quad = [3trp^-(trp3)2]/2. (9) 

-Pquad frequently used in the analysis of quantum mechanical density matrices. 
Another important purity measure for quantum density matrices is the von Neumann 
entropy, whose analog here is PvN, where 

1 ^ 

This is the natural purity measure in quantum statistical mechanics,^' and is con- 
sidered here for comparison. Each of the measures defined in Eqs. (IH]), ([9]), (fTOj) is a 
purity measure, and in a neighborhood of the totally unpolarized state, is linear in 
the eigenvalues. 

An important question is how the various 3-dimensional degrees of polarization 
measure a paraxially polarized beam; this is the case usually considered in discussions 
of Rayleigh scattering. In this case, the choice of cartesian coordinates in Eq. ([7]) 
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Fig. 2. Contours of constant P^y (thick lines) and Pquad (thin green lines) in 
the S*!, ^a-plane, for 3x3 paraxial partial polarization. The contour lines of 
-Pquad depend only on the radius >y Sf + 5*1, whereas the -pRay contours have 
weak angular dependence. 




(a) P2 (b) 

Fig. 3. Illustrating the various 3-dimensional polarization degree measures P3 
for paraxial partial polarization, parametrized by the paraxial measure P2- (a) 
P3 plot; (b) Plot of P3 — (P2 + l)/2 (i.e. linear part subtracted). P3 = Paay 
with 5*3 (Si) = (thick (dashed) line); Pun (thick gray line); Pquad (thin green 
line), PvN (thin dashed purple hne). 
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eliminates S2, so depends on the two Stokes parameters Si and S'3, namely 



1 

P3 ,trans o 



/ \ 

1 + ^1 -iSs 
iSs 1-Si 



y y 



(11) 



The transverse, two-dimensional degree of polarization is P2 = \/ Sf + (agreeing 
with the quantity defined in Ref. 4). However, as mentioned above, the transversally 
unpolarized state with Si = S3 = 0, has some measure of 3-dimensional polarization, 
as it appears completely unpolarized only in the propagation direction; its purity can 
neither be nor 1. P^ay is plotted as a function of 5*1 and 5*3 in Fig. [21 although its 
main dependence is on P2, there is weak dependence on the angle arctan(S'3/S'i). The 
paraxial partially polarized is pure when P2 = 1, and all definitions of polarization 



agree here. When P2 = 0, PRay = ^quad = ^im = 1/2, and P.^ = ^1 - log 2/ log 3 ^ 
0.608. The 3-dimensional degrees of polarization for this example are plotted in Fig. 

m 

Two different 3-dimensional polarization matrices with the same eigenvalues were 
considered in Ref. 3. Generalizations of those examples are as follows. The first is the 
diagonal matrix pg = diag{l + a, 1, 1 — a}/3, with < a < 1; this is completely 
unpolarized when a = 0. It is represented geometrically by an ellipsoid whose axes 
are the diagonal elements, and the mean angular momentum vector N is zero. The 
second example is defined such that each Mj = 1/3 (the ellipsoid is a sphere) and 
\N\ = N3 = a/3 {0 < a < 1). For example 2, 1(6, (f)) is uniform, and P2(^, 0)//(^, 0) = 
a\ cos6'|. The two examples have the same eigenvalues (1 + a)/3, 1/3, (1 — a)/3, but 
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geometrically are very different. In both cases, PRay depends linearly on a, with 

PRay,exl = 0.5932a, PRay,ex2 = 0.5a, (12) 

where the gradient for PRay,exi is determined by numerical integration. The two en- 
sembles thus have different Rayleigh-defined degrees of polarization for all a > 0. By 
comparison, unitary invariance means the three other measures do not distinguish be- 
tween the two ensembles, with Pun — 2a/3, and Pquad = a/y/S (Pvn is approximately 
linear in a, but cannot be expressed simply). 

The final example considered here is a geometrically simple case in which the 
polarization matrix is real and degenerate (its ellipsoid is axisymmetric and iV = 0) , 
i.e. pg = M = diag{2m, 1 — m, 1 — m}/2 for < m < 1. In this case, PRay can be 
found analytically, and PRay = Pun — Pquad — |3m — l|/2, and Pvn is numerically 
close. 

4. Discussion 

Elements of the polarization matrices are complex when the polarization states 

in the underlying ensembles have some elliptical or circular polarization; in pg, this 
is given by a nonzero value of the third Stokes parameter Ss, and in pg, by a nonzero 
net angular momentum vector N. Therefore, unitary invariance implies that states 
of elliptic polarization are treated equivalently to states of linear polarization. This 
is clearly appropriate for paraxially propagating states of polarization, where all po- 
larization states are on the same footing on the Poincare sphere, and may be freely 
transformed to each other by retarders and rotators without changing the degree of 
polarization. 
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The situation is less clear in nonparaxial physics, where there are several well- 
known deviations from the paraxial case, such as the polarization topology of in- 
homogeneous 3-dimensional fields, ^^'^^ and the geometric phase in twisted optical 
fibers. ^^'^^ This difference originates in the fact that the 3-dimensional orientation 
of a state of linear polarization requires two parameters (the direction of the real E 
vector), but the orientation of elliptic polarization requires three (the Euler angles 
of the polarization ellipse). In Rayleigh scattering, a ray's polarization is given by 
the projection of the 3-dimensional polarization matrix into the transverse plane of 
the ray by Eq. ([3]); complex states of elliptic polarization project over all directions 
differently than real states of linear polarization, and this is manifest in the lack of 
unitary invariance in P^ay, ultimately due to the non-unitarity of the physical mecha- 
nism of Rayleigh scattering. Given the geometric difference between nonparaxial and 
paraxial polarization, it seems that appealing to analogy with the two-dimensional 
paraxial case is insufficient motivation by itself for demanding the unitary invariance 
of the 3-dimensional degree of polarization; consideration of other physical instances 
of the 3-dimensional polarization matrix, such as in near-field diffraction, or light in 
a twisted fiber, may lead to other physically-motivated definitions of the polarization 
matrix that are, or are not, unitarily invariant. 

The definition can be modified mathematically to give a unitarily invariant 
version of pRay, at the cost of the physical interpretation. Instead of averaging over 
symmetric projections perpendicular to the observation direction u{6,(j)), one inte- 
grates over complex hermitian projections that are orthogonal to all complex vectors 
(i.e. all polarization states). This average respects unitary invariance, and can be 
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diagonalized before the average is taken. Therefore, under this unitary invariant mea- 
sure, the partially polarized paraxial ensembles behave like -PRay with 5*3 = (plotted 
in Fig. [3]), the examples 1 and 2 both behave like example 1, and any degenerate P3 
as the axisymmetric example. 

The geometrical interpretation of the 3-dimensional polarization matrix was used 
in Ref. 20 to describe the polarization of skylight using a canopy atmosphere model. 
Incident sunlight is Rayleigh scattered in a thin layer of fixed height above the plane 
of the observer, giving rise to P3 independent of position; the polarization in an 
observation direction was found using an orthogonal projection of this matrix (by Eq. 
(IH])); in particular, the neutral points, ^^'^^ where the degree of polarization vanishes, 
correspond to the Maxwell axes of the ellipsoid. 

Clearly, the definition of Paay in & is not the only way P2{9,(j)) and I{9,(j)) 
can be combined to give a well-defined degree of polarization; for instance, any 
positive powers can be used, and I{9, (p) need not be included. In particular, the 
choice P2{6, (j)Yl{9, 0)^ in the integrand can be easily integrated, yielding 2[3trM^ — 
(trM)2 + hN ■ A/']/[trM2 + 3(trM)2]. The choice of (ED is made here not only be- 
cause its definition is physically natural, but also it is closer, at least in the examples 
considered, to the measures and Pquad- In quantum mechanics, it is more usual 
to consider the squares of the purity measures considered in optics (such as the mea- 
sures in Eqs. ([9]) and ffTOl) ): the squares are more natural from the point of view of 
mathematical definition, but do not appear directly (in two dimensions) in terms of 
the elegant decomposition p2 = P2Ppure + P2p2,\in- 

It may be possible to determine the Rayleigh degree of polarization directly by ex- 
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perimental measurement. A possible Rayleigh scatterer in such an experiment would 
be a metallic nanoparticle, embedded in some transparent medium; this would be 
placed in an appropriate incoherent light field represented by the desired P3, suffi- 
ciently constant over the lengthscale of the particle. It would only be necessary to 
measure the polarization of p2{0,(f)) over a hemisphere, as the projections in Eq. ([3]) 
do not distinguish ±u{6, 0) - measuring the polarization over a large number of scat- 
tering directions may cause technical problems, and further discussion is outside the 
scope of the present work. Of course, the average over all directions of the scattered 
ray polarization depends on the kind of scatterer (anisotropic Rayleigh, Mie, etc); it 
is likely that different, scatterer-dependent measures of the degree of polarization can 
be defined. However, these other scatterers have some internal structure, unlike an 
isotropic Rayleigh scatterer, and the Rayleigh degree of polarization depends only on 
the electromagnetic polarization matrix p^. 

Projecting a density matrix on subspaces is reminiscent of Gleason's theorem^^'^^ 
in quantum mechanics, in which projection assigns weights to 1-dimensional sub- 
spaces. However, averaging such one-dimensional projections (over all hermitian pro- 
jections, being quantum mechanical) only gives the trace; this construction completely 
misses the complicated structure of pRay, and even its unitary counterpart. 

The 3-dimensional polarization matrix is insufficient to describe propagation of 
partially coherent nonparaxial fields, since propagation direction information is also 
required (in the form, for instance, of a Wigner function^^). It is not clear physically 
how a 3-dimensional degree of polarization ought to behave under propagation, espe- 
cially if it is assumed to depend only on the polarization matrix. Mechanisms such as 
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dipole scattering, however, are completely independent of the propagation directions 
of the electromagnetic fields giving rise to the polarization matrix (neglecting recoil 
effects), and are therefore a natural place to find physical relevance of measures of 
the 3-dimensional degree of polarization. However, as has been shown, the natural 
measure based on this process does not depend only on the eigenvalues of (it is 
not a unitary invariant), unlike other measures which have been proposed. The as- 
sumptions behind the definition of the 3-dimensional degree of polarization therefore 
require further consideration. 
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